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Abstract 



In the center of the Milky Way, as well as in many other galaxies, a compact star cluster 
(~| around a very massive black hole is observed. One of the possible explanations for the 

Qh formation of such Nuclear Star Clusters is based on the 'merging' of globular clusters in the 

Q inner galactic potential well. By mean of sophisticated N-body simulations, we checked the 

?H validity of this hypothesis and found that it may actually has been the one leading to the 

!/3 formation of the Milky Way Nuclear Star Cluster. 



1 Introduction 



> 

^" The existence of very compact star clusters around the centers of many galaxies is well 

^ ascertained, nowadays. These clusters, called Nuclear Star Clusters (NSCs) are among the 

C^ densest star clusters observed, with effective radii of a few pc and central luminosities up 

t — to 10 L0. Recently, the study of dense stellar systems has raised new interest because of 

^^ the discovery, in an ever increasing number of galaxies, of compact nuclei in form of stellar 

I resolved systems and/or NSCs which are now known to be present in galaxies across the 



JV whole Hubble sequence and not only in the dE, N galaxies (Bekki & Graham 2010). A 

• ^H radial profile and velocity structure for a NSC can be determined only for the Milky Way 

rS NSC which is close enough that it can be resolved into individual stars. The Milky Way NSC 

C^ has an estimated mass of 10 M©, and it hosts a massive black hole whose mass, 4.3 x 10^ 

Mq, is uniquely well determined. 

The formation mechanism of nuclear star clusters is unknown. Two competing models are 
possible. In the gas model, a NSC can form from the gas that migrates to the center of the 
galaxy where then forms stars. A variety of scenarios have been proposed to account for the 
required fast radial inflow of gas into the galactic center, including the magneto-rotational 
instability in a differentially rotating gas disk, tidal compression in shallow density profiles 
or dynamical instabilities. Alternatively, in the merger model, massive clusters migrate to 



the center via dynamical friction and merge to form a dense nucleus (Tremaine et al. , 1975 



Capuzzo-Dolcetta, 1993). Observations of NSCs in dE galaxies suggest that the majority, but 
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not all, dE nuclei, could be the result of packing mass in form of orbitally decayed globular 
clusters (GCs) ( |Lotz et al. 2004). Numerical simulations have also shown that the basic 
properties of NSCs, including their shape, mass density profile, and mass-radius relation are 



reproduced in the merger model under a variety of explored conditions (Capuzzo-Dolcetta et 
aH |2008a|b] [Hartmann et al.| [20lT| ). 



2 The Milky Way nuclear cluster 

As we said, a NSC is present in the central region of the Milky Way. For obvious reasons, it is 



the only NSC which can be resolved into individual stars in spite of its high density ( Schodel 



et al. 2007); moreover, detailed kinematic studies exist that allows a precise determination 
of the mass of the central black hole, Mbh = 4.3 x 10® M©. Moreover, the relaxation time 
at Sgr A* influence radius is robustly estimated as 20 - 30 Gyr (Merritt, 2010) suggesting 



there has not been time for a Bahcall-Wolf cusp to rise. This is consistent with the almost 



flat density distribution of late type stars within 0.5 pc from the center (Buchholz et al. 



2009). As stated in the Introduction, at least two competing models are possible for the 



formation of the MW NSC. These mechanisms must account for the observed evidence of the 
MW NSC, i.e.: a mass of the NSC M^sc ~ lO'^ Mq, a density profile with a core of about 
0.5 pc and decreasing as r~^-* up to r ~ 30 pc. Out of 30 pc there is a large nuclear stellar 
and molecular disk of approximately same radius. 

2.1 Globular cluster merging and MW NSC formation 

We modeled the final evolution of a set of massive GCs which decayed in the inner region 
of the Milky Way due to the dynamical friction braking exerted by the galactic stellar field. 
The details of the computations and the extended presentation and discussion of results are 
given in 



Antonini et al. (2011); here we give just a brief summary. 



Our Galaxy, likely triaxial in its inner part, is dense enough to make the orbits of GCs with 
mass around 10® Mq to decay (by dynamical friction) to the inner 50 pc in a time significantly 
shorter than a Hubble time. We thus consider 12 such GCs and their final evolution toward a 
merger state studied via N-body simulations of their motion in the inner region of the Milky 
Way, represented self-consistently as a set of N particles. The N body galactic environment is 
sampled from a shallow (p oc r^^") cusp; the contribution of the central massive black hole 
is also considered. Simulations are done with both the high precision, parallel PhiGRape 



(Harfst et al. , 2007) and NBSymple (Capuzzo-Dolcetta et al. , 2011) codes running on the 



RIT Grape cluster and on the CPU+NVIDIA TESLA C2050 platform in Roma, Sapienza. 



3 Results 



The merging occurs rather quickly: after about 20 crossing times the resulting system at- 
tains a quasi equilibrium configuration, as it is shown by the almost steady time profile of 
Lagrangian radii. This corresponds to a slowly evolving (via relaxation) super stellar cluster 
(SSC) which oscillates around the massive black hole. Figure [T] shows the growth with time 
of the spatial density profile after the various merger events. The core structure is conserved. 
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Figure 1: Spatial profile of the central NSC after 3, 6, 9 and 12 mergers. The central density 
grows with time. The dashed line is the fit to the NSC profile obtained at the end of the 
entire simulation using a broken power law model. 

Figure^ (left panel) gives the comparison between the surface density of the 12 GCs merger 
and the underlying galactic stellar profile. The kinematic properties of the merger system 
are resumed by the Fig. [2] (right panel), which shows the presence of a tangentially biased 
zone, intermediate between an inner (close to the massive BH, r < 0.3 pc) and outer (r > 40 
pc) isotropic regions. This tangential anisotropy grows with the number of infalls. This may 
be an important signature of the modes of NSC formation. 

The SSC has a density which is not far from the sum of the individual GC densities, 
so that the projected density profile in the nuclear region is remarkably similar to that of the 
Milky Way NSC. The final NSC morphology is nearly oblate with a value T ~ 0.2 for the 
triaxiality parameter. 

4 Conclusions 



As alrea dy shown by previous papers (Capuzzo-Dolcetta, 1993; Capuzzo-Dolcetta et al. 
2008a|b ) the dynamical friction dragging of massive globular clusters toward the galactic 
center is a viable explanation of the local high densities. The novelty of this work is the 
inclusion in the simulations of a massive black hole in the galactic center. Its role, although 
important in the vicinity of its influence radius, does not alter the general characteristics of the 
merger event of a set of GCs falling to the center. The merger reaches a sort of quasi-steady 
state, slowly evolving due to internal relaxation which is partly affected by the presence of 
the massive black hole. With regard to the Milky Way nuclear cluster, it is relevant noting 
how the original 'core' profile of the GCs (the NSC building blocks) is maintained for a time 
long enough to justify the core actually observed in the MW NSC. Moreover, the outward 
profile radial slope, p oc r~^, of the merger system is in good agreement (within 30 pc from 
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Figure 2: Left panel: projected density profile of the merged system at the end of the 
simulation (solid line). Dashed line refers to the galactic density profile. Right panel: radial 
behaviour of the anysotropy parameter /3 = 1 — cr|/a^. 

the center) with the observed MW NSC profile. Surely, an even more stringent answer to the 
question of the MW NC formation will come from a high resolution, well extended in time 
numerical study of the secular evolution of the merger remnant as well as from a detailed 
study of the expected observed stellar population (both in progress). 
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